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Abstract 

Additional experimental evidence is presented in support of the recent hypothesis that a possible solar influence could explain 
fluctuations observed in the measured decay rates of some isotopes. These data were obtained during routine weekly calibrations of 
an instrument used for radiological safety at The Ohio State University Research Reactor using 36 C1. The detector system used was 
based on a Geiger-Miiller gas detector, which is a robust detector system with very low susceptibility to environmental changes. 
A clear annual variation is evident in the data, with a maximum relative count rate observed in January /February, and a minimum 
relative count rate observed in July/August, for seven successive years from July 2005 to June 201 1. This annual variation is not 
likely to have arisen from changes in the detector surroundings, as we show here. 
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1. Introduction 

Evidence for a possible solar influence on nuclear decay rates 
has recently been presented based on the analysis of decay rate 
measurements taken at three independent institutions. The first 
was an apparent change in the measured decay rate of 54 Mn dur- 
ing a series of solar flares in December of 2006 [1]. The 54 Mn 
data were being collected as part of a half-life measurement 
utilizing continuous four-hour measurements. This allowed a 
time resolution capable of seeing changes that could have been 
caused by a solar flare, which typically lasts minutes to hours. 
This work was then followed by two additional papers by our 
group , where data were analyzed from half-life measure- 
ments taken by two independent groups, one at the Brookhaven 
National Laboratory (BNL) in Upton, New York, USA, and 
the other at the Physikalisch-Technische Bundesanstalt (PTB) 
in Braunschweig, Germany. The BNL group had undertaken a 
measurement of the half-life of 32 Si, and the data from that ex- 
periment exhibited a periodic oscillation with an approximate 
period of 1 year The measurements taken at the PTB in 
Germany were of a 226 Ra standard used for comparison in the 
measurements of the half-life of I52 Eu |5|, and in a longer-term 
analysis of the stability of detectors used in standards laborato- 
ries. The 226 Ra data also showed a periodic oscillation, again 
with a period of approximately 1 year. 

In the subsequent analysis of the raw data obtained from the 
BNL and PTB experiments, both data sets were shown to have 
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not only the same period, but in the two years during which the 
data sets overlapped they had the same approximate phase and 
amplitude as well 12j,|6J]. Moreover, it was shown that both data 
sets were not only in phase with each other, but also appeared to 
be approximately in phase with the distance of the Earth from 
the Sun. Taking all of these experiments into account (BNL, 
PTB and Purdue), a reasonable case could be constructed for 
the possibility of a solar influence on nuclear decays yfl. This 
case was subsequently strengthened as a result of an analysis 
by our group [7, 8], where an additional periodicity was iden- 
tified in the BNL data at 1 1.25+0.07/yr, and in the PTB data 
at 1 1.2 1+0.1 3/yr. Both of these peaks may be linked to the ro- 
tation (and probable inhomogeneous nature) of the core of the 
Sun |0- IT ] . A third periodicity was also identified by our group 
1 12] in both the BNL and PTB data sets, which is analogous to 
the Rieger periodicity lfl3ll with a period of approximately 173 
days. An analysis of the phases of these periodicities was also 
carried out IU411 which determined that the phase characteristics 
of the annual periodicities could reasonably be attributed to a 
solar influence on the decay rates. 

The suggestion of a solar influence on nuclear decay rates 
has been met with some criticism, however. An analysis by 
Norman et al. [ 15] of decay data taken for several isotopes in 
their laboratory did not see evidence for an annual effect similar 
to that reported in Ref. J2]. Silverman and Strange lfl6Tl exam- 
ined data from 22 Na decay measurements, and also found no 
similar fluctuations. Cooper [ 17] analyzed the heat output data 
from the radioisotope thermal generators (RTGs) onboard the 
Cassini spacecraft, and found no evidence for time variation in 
the decay of 238 Pu. Based upon the absence of oscillations sim- 
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ilar to those found in the BNL and PTB datasets, both Norman, 
et al., and Cooper concluded that the solar influence suggested 
in Refs. lUrjlt] was not present. However, a more recent analysis 
1 18|] of the data of Norman et al. B15I1 does suggest the presence 
of a solar influence, albeit at a lower level than indicated by the 
BNL or PTB data. As has been noted previously [19], the very 
same nuclear physics considerations which are responsible for 
the fact that beta-decay half-lives vary from fractions of a sec- 
ond to billions of years (i.e., nuclear wavefunctions, selection 
rules, phase space, etc.) would also apply to the effects of any 
solar influence. Hence, there should be no expectation that pe- 
riodic effects would be present in all beta-decays at the same 
~ 3 x 1(T 3 level seen in the BNL and PTB data. This observa- 
tion applies as well to the data in Ref. lfl6Tl . and especially to 
the analysis in Ref. [ 17]. In the latter case, the decay of 238 Pu is 
a pure alpha-decay, and leads to a daughter ( 234 U) which is an 
alpha decay with a 246,000 year half-life ll20Tl . Therefore, there 
would have been no significant contributions from beta-decays 
in the Cassini RTG data, whereas all previous data sets in which 
periodic effects were seen were measurements of beta-decays or 
electron capture. 

Similar anomalous behaviors and periodicities in nuclear de- 
cay data have in fact been observed by other groups. Data from 
the measurement of 60 Co and 9() Sr/ 90 Y published by Parkhomov 
1I21I - I23I1 also exhibit annual and monthly periodicities when 
measured on separate Geiger-Miiller (G-M) counting systems 
in a controlled experiment. Interestingly, 239 Pu counted by 
Parkhomov did not show any such periodicities. This demon- 
strates two important points: first, that the oscillations were 
likely not of an environmental origin; and second, the oscilla- 
tions appear to arise primarily in beta-decays, in agreement with 
the previous remarks. Also, since the counting systems were G- 
M detectors in the Parkhomov experiments, there would be no 
environmentally induced gain shifting since there is no ampli- 
fier in the system. 

Further evidence of annual periodicities in decay data was 
presented in an earlier publication by Ellis J3- His data 
showed an annual oscillation in the measured decay rate of 
neutron-activated manganese foils used to calibrate a system of 
plutonium-beryllium neutron sources. Interestingly, the 56 Mn 
counts exhibited an annual periodicity, yet the 137 Cs standard 
used to calibrate his scintillation detection system did not. This 
indicates that, as an experimental observation, isotopes have 
different sensitivities to whatever influence is causing the ob- 
served effects. Moreover, since the two isotopes were measured 
on the same counting system, it would also appear to rule out a 
simple environmental systematic cause. This supports the anal- 
ysis by Jenkins, Mundy and Fischbach [19], who examined all 
of the likely environmental influences on the counting systems 
utilized in the BNL, PTB and Purdue experiments J2.13.l5l] and 
concluded that all of the known suspect effects were too small 
to have caused the observed oscillatory behavior. Additionally, 
a recent paper by Steinitz et al. I25I1 presents extensive evidence 
for annual and sub-annual periodicities in the measured decay 
rates of 222 Rn (and its progeny), which lends further support 
to the solar influence hypothesis. Furthermore, we performed 
additional analyses j26ll of the 90 Sr/ 9l) Y data published in Refs. 



I21I-I23I1. and found that the data which contained annual and 
monthly periodicities also exhibited striking similarities to the 
frequency content in the Mount Wilson Solar Observatory's so- 
lar diameter measurement data. 

Some recent experiments have been performed to test the hy- 
pothesis of a neutrino-mediated solar influence on terrestrial 
nuclear decays by two independent groups, two of them per- 
formed by our group and one performed by an independent 
group in South Africa. The two experiments at the National 
Institute of Standards and Technology (NIST) performed by 
our group H27L 12811 examined a variation of the electron anti- 
neutrino flux (v e ) resulting from the /T-decay of 198 Au by uti- 
lizing samples with different geometries. The theory behind this 
experiment is that if the specific activity was high enough from 
the resulting decay of 198 Au in a sphere (compared to a foil as in 
Ref. [27] or a wire as in Ref. [28]), the v e flux could approach 
the solar electron neutrino (v e ) flux experienced on Earth, which 
is ~ 60 x 10 9 v e cmT 2 s~ x . Although the results of both of these 
experiments were inconclusive, they did not clearly support the 
hypothesis that v e could affect the 198 Au /T-decay. A differ- 



ent experiment was performed by de Meijer et al. 02911 . where 
decay rates of various isotopes were measured in close proxim- 
ity to a 2 MW,i nuclear reactor, which is a well characterized 
sources of v e . The negative results of this series of experiments 
led de Meijer et al. B29I1 to suggest two possibilities. The first 
was that the v e flux was not high enough from the 2 MW^ re- 
actor, and a larger reactor may show better results. The second, 
which was similar to the conclusions drawn by our group in 



Refs. 1 27, 28|> was that v e might not have the same effect as 
v e on f3~ -decay. Furthermore, the possibility exists that solar v e 
are not primarily responsible for the observed effect, but rather 
some other component of solar neutrino flux (e.g., v M or v T ) or 



an as yet unknown particle or field [28]. 

The purpose of this article is to present 36 C1 decay data col- 
lected at The Ohio State University Research Reactor (OS- 
URR), in Columbus, Ohio, USA, over the course of 7 years, 
which further strengthen the case for a solar influence on some 
nuclear decays. The data were taken weekly, as part of the cal- 
ibration check of an instrument used at the OSURR facility, 
thus the data were not the result of an experiment per se, but 
were collected as part of routine operations at the OSURR. It 
is evident from the data shown in Fig.[T]that there is an oscilla- 
tion with an approximate annual period that appears to correlate 
with the inverse-square of the Earth-Sun distance, and possibly 
additional frequencies that can also be linked to the Sun. 

2. Data Collection/Detector Set-up 

As noted above, the 36 C1 data were collected weekly as part 
of the efficiency check of the Eberline Beta Counter BC-4 that 
is used for counting of contamination survey wipes. The de- 
tector system incorporates a 1 .75 inch (4.4 cm) diameter G-M 
pancake-style tube contained inside 0.875 inches (2.22 cm) of 
lead to reduce background lf30tl . The BC-4 instrument itself sits 
in the reactor bay at OSURR, which has partial environmen- 
tal control. Although it is not air-conditioned, space heaters are 
used to maintain the interior temperature in a comfortable range 
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Figure 1 : Plot of the measured 36 C1 decays taken at the Ohio State University Research Reactor (OSURR). The crosses are the individual data points, and the blue 
line is an 1 1-point rolling average centered on each point. The red curve is the inverse of the square of the Earth-Sun distance. Error bars are shown for a limited 
number of 36 C1 points in order to maintain readability. 



during colder months. Thus, temperature is controlled to some 
extent, but relative humidity is not. 

The source utilized for the calibration check is a 0.4 /iCi 36 C1 
split check source (manufactured by Nuclear-Chicago, Model 
SK2-1) with a diameter of 1.0 inch (2.54 cm). The two alu- 
minum half-disk sources which comprise the check source are 
contained within an aluminum holder with an outer diame- 
ter of 1.25 inches (3.2 cm). The active regions of the source 
are two machined circular depressions, each approximately 0.2 
inches in diameter, located near the center of the 1-inch disk on 
each half disk. The decay of 36 C1 is primarily by p~ -emission 
(98.1%, £,,=708.6(3) keV, T l/2 = 3.01(2) x 10 5 y HI) to the 
ground state of 36 Ar, which is stable (there is also a competing 
K-capture mode to 36 S with a 1.90% branching ratio). We note 
in passing that this is the same isotope used as the standard in 
the BNL 32 Si half-life experiment. 

After early 2006, the source geometry within the BC-4 
counter was controlled by an aluminum disk insert with an outer 
diameter of approximately 2.0 inches. This disk centered the 
source within the 2-inch planchet used to place the source on the 
detector tray, and this ensured that the check source was located 
in the same position under the detector for each calibration- 
check measurement. The source/insert combination was then 
placed into the planchet, and then the planchet was placed into 
the source tray in the BC-4, which slides in to locate the source 
under the detector within the counter. Geometry in the counter 
is well controlled, such that the source-detector separation is 
very small (~2 mm). The slide tray position is controlled by 
a stop at the rear of the detector system, thus controlling the 
y-axis geometry. 



Data recorded each week were the gross counts (source with 
background, with no dead-time correction), and a ten minute 
background count. Typical 36 C1 count rates were approximately 
124,600/120s (~ 1038 counts/s), with a fractional uncertainty of 
~0.28%. Typical background count rates were approximately 
225 counts/600s (~0.35 counts/s), with a fractional uncertainty 
of ~6.7%. Since the half-life of 36 C1 is much longer than the 
duration of the data series reported here, the standard decay rate 
can be assumed to have been approximately constant during the 
course of the experiment. 



3. Results/Discussion 

The measured 36 C1 data are shown in Fig. Q] These data 
points are the weekly counts from 7 January 2005 to 17 June 
201 1, a total of 334 points. Since the 36 C1 decay rate can be as- 
sumed to be relatively constant over the five years of the series 
reported here, the only adjustment to the data presented in Fig.Q] 
is to normalize each measured count by dividing the measured 
counts by the average of all 334 counts (x = 124593, <x = 632). 
Note that the consistency of the data improves after early 2006 
when the source geometry control was added. An 1 1 point mov- 
ing average of the data is has been added to the figure as a solid 
blue curve to aid in the visual presentation. 

It is clear from Fig.Q]that the data exhibit an annual periodic- 
ity, with higher counts in winter (approximately February), and 
lower counts in the summer (approximately August). This pe- 
riodicity is very similar to the annual oscillations presented in 
the 32 Si/ 36 Cl data from BNL HSl and the 226 Ra and progeny 
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data from the PTB [2, 5]. Further analysis reveals that the nor- 
malized counts shown in Fig. Q] oscillate with an approximate 
annual period that is roughly correlated with the annual period 
of the inverse-square of the Earth-Sun distance, which is shown 
as the red curve in Fig.Q] 

To verify the presence of the annual periodicity, a power- 
spectrum analysis was performed for the data in Fig. Q] In 
power-spectrum analysis, the probability of the null hypothe- 
sis (that the peak is due to normally distributed random noise) 
is given by e~ s (where S is the power value of the frequency 
peak) 13211 . These results are presented in Fig. |2l with fre- 
quency in units of inverse years. The prominent peak in the 
figure for a frequency of 1 yr _I (with a power of ~54) con- 
firms the annual periodicity. The peak in the power spectrum at 
~0.25 yr" 1 is probably related to the length of the data set, but 
may represent the influence of the ~1 1-year solar cycle. None 
of the other peaks in the power spectrum are statistically signif- 
icant. Although the annual periodicity shown in Fig. [2]is quite 
prominent, there remains a question as to the cause of this an- 
nual oscillation. Is the annual periodicity related to some prop- 
erty of the Earth-Sun system, such as the inverse-square of the 
Earth-Sun distance, or rather is it a consequence of interfering 
or modifying inputs to the counting system that vary periodi- 
cally with a frequency of 1 yr -1 ? Background radiation is an 
example of an interfering input, and the dependence of the de- 
tector, counting system, and the transport of radiation across the 
source/detector gap on environmental factors (such as tempera- 
ture and humidity) is an example of a modifying input. 
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Figure 2: Frequency spectrum for the OSURR 36 C1 data. The annual oscillation 
is clearly evident, as is a peak at ~0.25/yr that is probably related to the length 
of the data set. None of the other peaks are statistically significant. 

To address these concerns, we begin by analyzing the de- 
tector system utilized in this experiment, which is G-M tube as 
described previously. A G-M tube is a gas detector with an inert 
gas (usually argon or neon) quenched with a halogen gas in this 
case. The use of a halogen quench over an organic compound 
as the quenching gas allows for a recombination of the diatomic 
halogen, and extends the life of the tube indefinitely. When ra- 



diation, such as an a- or yS-particle, or an energetic electron 
from a photon interaction in the detector wall, interacts within 
the tube gas, an electron-ion pair is generated. Photons can in- 
teract directly with the working gas and cause ionization, but 
this is a lower probability interaction. The bias of the detector 
provides a driving force for the electrons to move to the anode, 
and the ions to move to the cathode. As an electron is acceler- 
ated by the very strong electric field within the detector, it can 
accelerate to the point where it generates additional ionizations, 
which result in an avalanche of ionizations within the working 
gas. The avalanche of ionizations continues to the point where 
it self-terminates due to a reduction of the electric field due to 
the presence of so many positively charged ions. As a result, es- 
sentially all of the gas that can be ionized is ionized with each 
discharge. 

The large number of electrons collected at each discharge 
provide a charge pulse of sufficient magnitude that it does not 
require amplification to be counted, and therefore there are no 
pre-amplifiers or amplifiers required in the processing of the 
signal out of the detector. A trigger circuit is the next step in 
the system, and as long as the pulse out of the detector is of 
large enough magnitude, a standard logic pulse is output, which 
is counted by a sealer-timer. While the disadvantages of the G- 
M detector are mostly along the lines of its inability to perform 
energy-spectroscopy or discrimination between radiation types, 
the significant advantage is that it is a robust, very stable sys- 
tem for detecting radiation events. Therefore it is suitable for 
relatively harsh working environments, and is a very practical, 
stable system. The general stability of G-M systems, similar to 
the one used to collect these data, with respect to temperature 
has been well studied II33I439I1 . These analyses determined that 
the stability of the G-M system was not likely to be a factor 
within the range of temperatures that were likely to be expe- 
rienced in the OSURR laboratory (65-85°F, though the actual 
temperature range experienced was likely much smaller). We 
can thus conclude that a change in the detector system is not a 
likely candidate for the source of the fluctuations. 

In light of the fact that beta-particles lose energy with each 
interaction with an air molecule, the transport of the beta- 
particles across the source-detector gap clearly depends on the 
density of the intervening air. This very argument, that air den- 
sity was a factor in changing the measured count rates, was one 
of the bases proposed by Semkow et al. B40I1 to question the 
claims in Ref. [2]. Semkow et al. J40ll assumed that humid 
air, such as would be experienced in the summer months, was 
denser than dry, winter air. This is not the case, as shown in 
Fig. [3] which is a plot of air density as a function of temper- 
ature and moisture content (measured by relative humidity) as 
adopted in CIPM-2007 liulial . As can be seen from the fig- 
ure, cooler, drier air is significantly denser than warm, humid 
air. Therefore, based on air density, if the count rates were to 
be affected by the energy loss of the beta particles due to denser 
air, one could assume that the count rates should be lower in the 
winter when the inside air is coolest and driest. Examining Fig- 
ureO] this is clearly not the case. In fact, the counts are actually 
highest when the air is the coolest and driest (i.e., the highest 
density). 
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Figure 3: Plot of air density as a function of temperature and relative humidity 
as adopted by CIPM |4l|.[4^1. 



It is still worthwhile to examine the beta particle transport 
across the source-detector gap using a radiation transport code, 
such as MCNP (something that was suggested as well in Ref. 
1 40]), even though the higher count rates in the cooler, drier, 
denser winter air imply that beta particle transport is not an 
issue across the very small gap between source and detector. 
Such an MCNP analysis was performed for the BNL source- 
detector system as a part of a larger, detailed analysis of envi- 
ronmental influences by our group on the BNL and PTB detec- 
tor systems in Ref. [19]. The BNL source/detector system (a 
gas proportional detector) was modeled using the Monte Carlo 
N-Particle extended (MCNPX, H) code for both the 36 C1 and 
32 Si samples. The 36 C1 source-detector distance in the BNL ex- 
periments was 4.000 mm, which is almost twice the distance 
between the source and detector for the OSURR BC-4 (2 mm). 
The results for the 36 C1 analysis from Ref. [19] are presented 
in Table [3] For each of the parameters analyzed, it can be 
seen from these results that any possible variation (measured 
in change per degree Fahrenheit, or A/°F) of parameters (e.g., 
energy deposition in the detector volume, particle current, etc.) 
is at least two orders of magnitude too small to have caused the 
observed periodic effects. This supports the conclusions of Al- 
burger et al. J3l in their own analysis of possible environmental 
and systematic influences on their detector system. In Table 1 
of Ref. |4J], limits were also set on other systematic influences 
at less than a 10~ 4 fractional change for each parameter, which 
is too small to have explained the oscillations they observed. 
We can thus safely say, based on the above analysis, that air 
density changes were not likely to have affected the transport 
of particles across the 2 mm source-detector gap, nor should 
it have affected the deposition of energy by the beta particles 
in the detector volume. This leaves one additional parameter 
of a gas detector system, which is the energy needed to ionize 
the working gas. Since it is known that a G-M detector only 
requires a single ionization to discharge essentially the entire 
tube, and also that the ionization potential of argon (the work- 
ing gas) does not change for pressures between 1 and 200 bars 
1 44], we can again safely conclude that a G-M tube is insensi- 
tive to small changes in the ambient laboratory pressures. 



Table 1 : Table of MCNPX results showing the sensitivities for various parame- 
ters of the BNL source-detector system, per degree Fahrenheit. [19?] 

Per Source e~~ Norm.(70°F) 
Det. E Deposition (MeV/ptcl) (A/°F) 0.34(40) X KT 6 39.8(67) X HP 

Det. Window e- Current (A/ F) 20.3(13) X lO" 6 35.3(33) X Ur 6 

Det. Window ECurrent(Me V/ptcl) (A/°F) 4.32(32) X lO" 6 32.4(34) X 10" 6 
Det./Source e~ Current Ratio (A/°F) 29.5(24) X lO" 6 41.5(48) X 10" 6 

Det./Source E Current Ratio (A/°F) 24.9(25) X lO" 6 33.7(48) X 10 -6 



Further analysis of climate as a cause of the observed oscilla- 
tions can be carried out by examining the phase of the 36 C1 data 
taken at OSURR, and comparing it to the phases of temperature, 
atmospheric pressure, and relative humidity data reported for 
Columbus, Ohio. The dominant frequency in the power spec- 
trum of 36 C1 provided in Figure [2] is at 1.06+0.09 yr _1 , and it 
has an associated phase of -48.8+3.2°, which corresponds to 
approximately 17 February, with 1 January equal to 0°. This is 
consistent with the results of frequency and phase analyses car- 
ried out by our group for the other experiments (BNL, PTB and 
that from Ref. [24]) presented in Refs. It is clear from 

Fig. Q] that the 36 C1 phase lags the phase of 1/R 2 , but this was 
also the case for the BNL [4] and PTB [5] examined in Ref. lyy, 
in which we showed that the phases of the BNL and PTB data 
were similar for the two year period in which the experiments 
overlapped. One can also see from Fig. Q] that the phase lag is 
not perfectly consistent across the six year period reported, a 
point to which we will return shortly. Thus, the reported phase 
of the 36 C1 data, at -48.8+3.2°, is an average, and gives no rep- 
resentation of phase stability. 

We next turn to the comparison of the 36 C1 data phase to the 
phases of the ambient temperature, pressure, and relative hu- 
midity in Columbus, Ohio, as reported by the National Climac- 
tic Data Center (NCDC) which is part of the National Oceanic 
and Atmospheric Administration (NOAA). The phases of am- 
bient temperature, atmospheric pressure, and relative humid- 
ity were -16.7°, -129°, and -143.6° respectively. We should be 
careful to note here that the phases we are reporting for the cli- 
mate data were obtained by sampling the climate parameters at 
the same sampling rate and dates as the 36 C1 data (i.e., the data 
reported for the days which correspond to the days of the BC-4 
calibrations). We should also note that the determined phases 
for the environmental variables are also averages, as for the 36 C1 
data, and do not give any indication of phase stability. In this 
case, as with the analysis by our group presented in Ref. |0] we 
see that the phases of the weather variables or their reciprocals 
do not appear to match the phase of the measured 36 C1 data. 

Although the phase of the 36 C1 data does not match that of 
1 /R 2 exactly, this question was addressed by our group in Ref. 
1 14], where the difference is attributed to an inhomogeneous 
solar core. This inhomogeneity then results in a non-uniform 
solar v e production between the northern and southern solar 
hemispheres. Earth's maximum exposure to the southern hemi- 
sphere of the Sun is ~8 March due to the tilt of the solar axis 
with respect to the plane of the Earth's orbit. Thus, the com- 
bination of the closer proximity at perihelion, and the differing 
exposures to more active regions of the Sun, could explain the 
phase lag. This phase lag is, in fact, very close to what is ob- 
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served in data from two of the major neutrino detectors. Anal- 
ysis of data from Super-Kamiokande-I in Japan |45, 46j] pro- 
vides a peak in the solar neutrino flux at perihelion plus 13+17 
days l45l,l46ll . and places a 68% confidence range for the actual 
peak that covers a range that includes the first week of Febru- 
ary. Analysis of data from the Sudbury Neutrino Observatory 
(SNO) |47, 48] shows the best fit for peak solar neutrino flux 
at perihelion plus 40 days (~12 February) B48I1 . which is within 
2<x (<r «3.25 days, see above) of the calculated phase of the 
OSURR data, or about 17 February. 

Moreover, as described in Refs. the actual observed 

phase, if any, in nuclear decay data will be dependent on the 
sensitivity of the decaying isotopes to the external effect, and 
whether or not the isotope(s) being measured are parents or 
daughters. Clearly, if the isotope being measured is a second 
or third generation progeny of a radioactive parent, any pertur- 
bations to the decay rate of the parent will have to work their 
way down the decay chain to the daughter, grand-daughter, etc. 
Evidence of this fact was reported in the analysis of the BNL 
data 0] in Refs. 10 H. The 36 C1 data and the 32 Si (in which 
case, the 32 P daughter was the nuclide measured) did not exhibit 
the same phases, even though they were measured in alternating 
counts, where each was counted for thirty minutes respectively, 
ten times during a ten hour run. Another example of differing 
isotopic sensitivities was noted by Ellis |24J], where the 56 Mn 
data collected on a sodium iodide detector exhibited an annual 
oscillation, but the 137 Cs calibration source measurements taken 
on the same detector on the same days did not. 

Based on the considerations in the preceding paragraphs, 
we have shown that temperature variations are not a likely 
explanation for the observed periodicity. However, to be 
thorough we can examine the possibility that any of the 
materials of the BC-4 housing, the source, or the source 
holder may have expanded or contracted due to thermal ex- 
pansion. The materials which comprise the source tray 
and detector are steel (which has a linear expansion co- 
efficient of 7.3xKT f> in-in- 1 -°F- 1 , BxlO^m-nr'-K- 1 ) and 
aluminum (which has a linear expansion coefficient of 
^JxlO^in-in- 1 - ?- 1 , 22.2xKr 6 m-nr 1 -Kr 1 ). Since the alu- 
minum thermal expansion coefficient is larger, we can use it as 
a bounding value and assume that any expansion in steel will 
be smaller. 

Examining the simplest conservative case, then, we assume 
that the steel tracks that hold the source tray expand such that 
the source-detector distance was changed. We will examine the 
system in three dimensions, with the xy-plane being the source 
holder and the z-axis running from the origin at the center of 
the source up through the detector. Examining a change of the 
z-axis first, and using the (larger) aluminum expansion coeffi- 
cient, we find that for an initial source-detector distance of 2 
mm and a 10°C (~18°F) change in temperature, we wouldreal- 
ize only a 0.0004 mm change in the source-detector distance. 

Based on this analysis, we find that the expected material ex- 
pansion in the z-direction could not lead to a count rate change 
similar to those observed, particularly since the expectation 
would be a linear dependence, and not an r 2 , since we have 
what amounts to a 2n geometry. Furthermore, we can also ex- 



amine the possible shifts that would occur in the xy-plane as a 
result of thermal expansion. If we assume similar changes to 
the tray width and length that we expect on the z-axis, we can 
conclude that the expected changes would be similarly incon- 
sequential. The result would be a shift on the order of ~0.0202 
mm of what amounts to a small active source region (~ 0.4 cm 2 
combined area) under what is, as we described previously, es- 
sentially a 2ji geometry of the pancake G-M detector, which has 
an active window of ~ 15.5 cm 2 . This too will not have enough 
impact to cause the observed variation in count rates. 

Having addressed the likely suspects for systematic and cli- 
mactic influences, there remains the question of counts due to 
background. While the detector itself is well shielded, there 
will still be background counts included in the 2 minute gross 
count measurements. As noted above, 10 minute background 
counts were taken each week along with the 36 C1 calibration 
count. The average of the weekly background counts for the 
period January 2006-June 2011 is 0.37+0.14 counts/s. When 
compared to the 1038+3 counts/sec for the 36 C1 source, it is 
evident that changes in the background were not likely to be 
responsible for the decay rate variations. A point by point 
comparison of the 36 C1 counts/sec to the measured background 
counts/sec is presented in Fig. [4] As can be seen from this fig- 
ure, while the annual periodicity is clearly evident in the 36 C1 
data, there is no similar periodicity evident in the background. 
This is verified by a power spectrum analysis similar to the one 
performed on the 36 C1 data, which is presented in Fig. [5] No 
similar frequency structure is evident in these background data. 

The measured background counts are dependent on what was 
happening in the laboratory area surrounding the detector at the 
time each count was taken, such as the reactor power at the 
time of the background count, and what radioactive materials 
happen to be located near the detector. For instance, the obvi- 
ous spike in the background counts during the summer of 2008 
was the result of the storage of the stainless steel reactor control 
rods (which are radioactive as a result of neutron-activation) 
nearby during a routine periodic reactor maintenance outage. 
Note from Fig. [4] that there was no discernible increase in the 
36 C1 count rate during that period; in fact, this increase in back- 
ground occurred in the middle of the summer decrease in the 
36 C1 counts. Providing additional support is the analysis of Ref. 
lfl9h . which includes a rigorous examination of background ef- 
fects in all types of gas detectors (specifically, the BNL and 
PTB detector systems). It was determined that in all cases, the 
dose rates to the detector from known backgrounds (such as 
cosmic radiation) were too low to have accounted for the os- 
cillations in the BNL or PTB data, and by extension we can 
assume the same applies here. 

There are some key conclusions that can be drawn from this 
analysis of these background data. First, changes in background 
could not have been the cause of the oscillations evident in 
the measured 36 C1 decay rate during the recording period. The 
background data do not contain the same periodicity as the 36 C1 
standard measurements as seen in Figs. [2] |4] and [5] Moreover, 
even a ~480% increase in the background count did not ap- 
pear to discernibly affect the measured 36 C1 decay rate, with the 
elevated background still only comprising ~0. 18% of the to- 
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Figure 4: Comparison of the measured CI and background measurements (in counts/second) on the OSURR detector from January 2006 to June 2011. Repre- 
sentative error bars shown are standard Poisson ( ^fN) uncertainties for the respective measurements. The large increase in the background countrate in 2008 was 
associated with the reactor control rods being temporarily stored near the detector system during a routine maintenance outage. 



tal counts (1034 counts/sec for the 36 C1, 1.83 cts/sec for back- 
ground). Second, the absence of periodicities in the background 
similar to those shown the 36 C1 data address two of our specific 
questions. One, the stability in the background data speak to 
the relative stability of the counting systems with respect to the 
influence of the environment on the detector or source. Two, 
the observed fluctuations are not the result of transients in the 
detector electronics, since no similar fluctuations appear in the 
background counts which were taken immediately before or af- 
ter the 36 C1 counts. We can now reasonably conclude that what- 
ever is causing the count rates is not likely to be in the detector, 
or the environment surrounding the detector. 

4. Conclusions 

Having eliminated all of the possible known systematic 
causes of the fluctuations, we are left with the possibility that 
the measured decay rate changes are caused by an influence not 
of terrestrial origin, such as the Sun as suggested in Refs. |l|- 
S EH EH- There is now apparent corroboration for the 
annual periodicity in the 36 C1 data from the BNL experiment 
in Refs. [2, 4]. While the sub-annual periodicities present in 
the BNL and PTB data sets as shown in Refs. 0, El EH do 



not appear in this OSURR data, this is not surprising based on 
the statistics of the relatively small count rates (fractional er- 
ror ~0.3%). In contrast, the 36 C1 data of the BNL experiment 
did show additional periodicities Jj, |7J, [12J], but the statistics 
of those counts were much better, with each data point con- 
sisting of ~1.2xl0 6 counts (fractional error ~0.03%). How- 




Frequency (year ) 



Figure 5: Frequency spectrum for the OSURR measured background data from 
January 2006 through June 201 1 . 



ever, the annual oscillation is clearly evident in the OSURR 
data, and thus lends support to the conclusions of Refs. 
[Hl3,H[l2.[l3.l25D, while providing additional evidence of the 
annual effect in 36 C1 from more recently measured decay data. 

There is an obvious need for more experiments examining a 
wide variety of detector technologies and a large array of iso- 
topes utilizing appropriate controls, including directly measur- 
ing temperature, pressure and humidity effects on detector sys- 
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tems in environmental chambers. There is also a need to ex- 
amine records of experiments and calibrations such as this one, 
where data are available that could be analyzed for the pres- 
ence of unexplained oscillations in what should be randomly 
distributed data. A concerted effort will be required by many 
to address the questions that remain about the causes of these 
fluctuations, including the development of a physical model. 
Hopefully this will also lead to more accurate determinations 
of half-lives for the ~3300 known radioactive isotopes, and a 
better understanding of the physical world. 
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